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Abstract. Having conducted a search for the A~1.3cm (22 GHz) water vapor line towards galaxies with nuclear 
activity, large nuclear column densities or high infrared luminosities, we present H 2 O spectra for NGC2273, 
UGC5101, and NGC3393 with isotropic luminosities of 7, 1500, and 400 Lq. The H 2 O maser in UGC5101 is by 
far the most luminous yet found in an ultraluminous infrared galaxy. NGC 3393 reveals the classic spectrum of 
a ‘disk maser’, represented by three distinct groups of Doppler components. As in all other known cases except 
NGC4258, the rotation velocity of the putative masing disk is well below 1000hms”* * ** ^. Based on the literature 
and archive data. X-ray absorbing column densities are compiled for the 64 galaxies with reported maser sources 
beyond the Magellanic Clouds. For NGC 2782 and NGC 5728, we present Chandra archive data that indicate the 
presence of an active galactic nucleus in both galaxies. Modeling the hard nuclear X-ray emission, NGC 2782 
is best fit by a high energy reflection spectrum with Nh^IO^'* cm“^. For NGC 5728, partial absorption with a 
power law spectrum indicates A^h^8x 10^® cm“^. The correlation between absorbing column and H 2 O emission 
is analyzed. There is a striking difference between kilo- and megamasers with megamasers being associated with 
higher column densities. All kilomasers (Lh 2O<10L©) except NGC 2273 and NGC 5194 are Compton-thin, i.e. 
their absorbing columns are <10^'*cm“^. Among the H 2 O megamasers, 50% arise from Compton-thick and 85% 
from heavily obscured (>10^®cm“^) active galactic nuclei. These values are not larger but consistent with those 
from samples of Seyfert 2 galaxies not selected on the basis of maser emission. The similarity in column densities 
can be explained by small deviations in position between maser spots and nuclear X-ray source and a high degree 
of dumpiness in the circumnuclear interstellar medium. 
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1. Introduction 

The JKaKc = 6 i 6 — 523 line of ortho-H20 at 22 GHz 
(A~1.3cm) is one of the strongest and most remarkable 
spectral features of the electromagnetic spectrum. The 
transition, connecting levels located approximately 645 K 
above the ground state, traces warm (Tkin^dOOK) and 
dense (n(H2)^10^ cm“^) molecular gas (e.g. Kylahs & 
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Norman 1987, 1991). Observed as a maser, the line can 
reach isotropic luminosities in excess of 10^ L© (Barvainis 
& Antonucci 2005). Being traditionally known to trace 
oxygen-rich red giant stars and to pinpoint sites of (mas¬ 
sive) star formation, it has more recently become a major 
tool to determine geometric distances and three dimen¬ 
sional motions of nearby galaxies and to elucidate the nu¬ 
clear environment of active galaxies, allowing us to map 
accretion disks and to determine masses of nuclear engines 
(for recent reviews, see Greenhill 2002, 2004; Maloney 
2002; Henkel & Braatz 2003; Morganti et al. 2004; Henkel 
et al. 2005a; Lo 2005; for 3-D motions, see Brunthaler et 
al. 2005). 

Extragalactic H 2 O masers have been observed for a 
quarter of a century (Ghurchwell et al. 1977) and emis¬ 
sion has been reported from more than 60 galaxies (e.g.. 
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Henkel et al. 2005b, Kondratko et al. 2006). The masers 
can be classified as related (1) to star formation, (2) to nu¬ 
clear accretion disks (‘disk-masers’), (3) to interactions of 
nuclear jet(s) with ambient molecular clouds or to ampli¬ 
fication of the jet’s seed photons by suitably located fore¬ 
ground clouds (‘jet-masers’) and (4) to nuclear outflows. 
H 2 O masers with apparent (isotropic) luminosities 
< 10 L 0 , often associated with sites of massive star for¬ 
mation, are referred to as ‘kilomasers’, while the stronger 
sources, known to arise from the innermost parsecs of their 
parent galaxy, are referred to as ‘megamasers’. The statis¬ 
tical properties of the sample of extragalactic H 2 O masers 
are not well explored. The first and last such comprehen¬ 
sive study, that of Braatz et al. (1997), dates back to when 
little more than a dozen sources were known. 

In order to improve our knowledge of the statistical 
properties of the greatly enlarged maser sample, we an¬ 
alyze the correlation between maser emission and X-ray 
absorbing column density toward the nuclear engine. In 
Sects. 2-3 we summarize recent observations and show 
spectra of sources detected in this work. In Sect. 4, we 
present X-ray column densities for the entire extragalac¬ 
tic maser sample, including a detailed analysis of Chandra 
archive data of the megamaser galaxies NGC 2782 and 
NGC 5728. Specific statistical properties of the sample in¬ 
volving X-ray absorbing columns and maser luminosities 
are discussed in Sect. 5. Sect. 6 summarizes the main re¬ 
sults. 

2. Observations 

Our measurements were obtained with the 100-m tele¬ 
scope of the MPIfR at Effelsberg in June 2001, October 
2004 and January and February 2005. The full beam width 
to half power was 40". A dual channel HEMT receiver pro¬ 
vided system temperatures of 100-350 K on a main beam 
brightness temperature scale. The observations were car¬ 
ried out in a dual beam switching mode with a beam throw 
of 2' and a switching frequency of ~1 Hz. The autocorre¬ 
lator backend was split into eight bands of width 40 or 
80 MHz and 512 or 256 channels. Each band could be 
shifted individually by up to ±250MHz (±3000kms“^) 
relative to the recessional velocity of the galaxy. Ghannel 
spacings were ~1 and ~4kms“^. 

The pointing accuracy was better than 10". 
Galibration was obtained by measuring NGG 7027, 
3G286, or W3(OH) (for flux densities, see Baars et al. 
1977; Ott et al. 1994). The data were reduced using the 
GLASS and GREG packages of the GILDAS software. 

3. Results 

Table □ summarizes observations of our sample of galax¬ 
ies showing nuclear activity, high nuclear column densi¬ 
ties (see Sect. 4) or high infrared luminosities. Among the 
18 sources observed, three were detected. Figs.^Elshow 
the line profiles of the detected targets. Table [3 provides 
line parameters obtained from Gaussian fits. Adopting a 


Hubble constant of Hq = 75kms“^ Mpc“^, to be consis¬ 
tent with previous maser studies and to be close to the 
parameters of the standard A-cosmology (e.g. Spergel et 
al. 2003), isotropic luminosities are LH 2 O 7, 1500, and 
400 Lq for NGC 2273, UGC5101 and NGC 3393, respec¬ 
tively. 



Fig. 1. The A~1.3cm (22 GHz) H 2 O megamaser profile 
toward UGC5101, with a channel spacing of 4.2kms“^. 
02000 = 09** 35“ 51 !6, J2000 = 6H 2T 11". Velocity scales 
are with respect to the Local Standard of Rest (LSR) 
and use the optical convention that is equivalent to cz. 
I4ys = czsys = 11809 km s“^ (NASA/IPAC Extragalactic 
Database (NED)). IAsr“Lhel = ±3.69kms“^. 


3.1. UGC5101 

UGC 5101 hosts one of the most luminous megamasers 
yet detected. Along with NGC 6240 (Hagiwara et al. 2002; 
Nakai et al. 2002; Braatz et al. 2003), UGC 5101 contains 
only the second H 2 O maser so far encountered in an ultra- 
luminous infrared galaxy (ULIRG). Interestingly, neither 
of the two ULIRGs with H 2 O maser emission is known to 
host an OH megamaser. 

As indicated in Sect. 1, all H 2 O megamasers studied 
with milliarcsecond resolution have been found to be as¬ 
sociated with the active galactic nucleus (AGN) of their 
parent galaxy. Is this also the case for ULIRGs, which are 
characterized by extremely high rates of massive star for¬ 
mation? While masers associated with star formation are 
not known to reach luminosities in excess of 10 Lq (for the 
highest luminosities reached, see Tarchi et al. 2002b), the 
superposition of a large number of such sources might lead 
to cumulative luminosities that rival those of luminous nu¬ 
clear megamasers. There is, however, evidence against this 
case. A large number of masers with luminosities of 0.01- 
ILq should lead to detectable H 2 O emission not only in 
UGC 5101, but also in Arp 220. Because of the large num¬ 
ber of such putative star forming H 2 O maser sources, their 
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Table 1. H 2 O maser observations 


Source 

R.A. Dec. 

(J2000) 

14ys 

(cz) 

Type®) 

p)) 

-^FIR 

(log) 

M—range 

(km s“^) 

rms 

(mjy) 

Channel 

Width 

(kms“^) 

Epoch^^ 

(lO- 

K’ 

'Am-") 

NGC 1667 

04 48 37.1 

-06 19 12 

4547 

Sy 2 

11.0 

4000- 5100 

10 

1.1 

0205 

> 10 

B; B99 

NGC2273 

06 50 08.7 

-h60 50 45 

1840 

Sy 2 

10.1 

1400- 2400 

9 

1.1 

0105 

> 18 

X; Gua05 

IRAS 07145-2914 

07 16 31.2 

-29 19 29 

1698 

Sy 2 

10.0 

1100- 2200 

77 

1.1 

1004 

> 100 

B; B99 







1050- 2200 

35 

1.1 

0105 



UGC5101 

09 35 51.6 

-\-61 21 11 

11809 

ULIRG 

11.9 

10550-12800 

2 

4.2 

1004 

13 + 2 

C; Arm04 

NGC 2992*’) 

09 45 42.0 

-14 19 35 

2311 

Sy 1.9 


1770- 2900 

11 

1.1 

0205 

~ 0.1 

B; Gua05 

NGC 3393 

10 48 23.4 

-25 09 43 

3750 

Sy 2 

10.4 

2880- 4800 

14 

4.2 

0105 

44+25 
— 11 

X; Gua05 

IRAS 11058-1131 

11 08 20.3 

-11 48 12 

16437 

Sy 2 

11.2 

15500-17500 

7 

4.2 

0105 

> 100 

B; RisOO 

NGC 4418 

12 26 54.6 

-00 52 39 

2179 

Sy 2 

10.9 

1550- 2900 

33 

1.1 

1004 

> 10 

C; Mai03 

Mrk231 

12 56 14.2 

-h56 52 25 

12642 

ULIRG 

12.2 

12050-13200 

7 

1.1 

0205 

~ 20 

X,B; Bra04 

NGC 4939®> 

13 04 14.4 

-10 20 23 

3110 

Sy 2 

10.3 

2500- 3770 

8 

1.1 

0105 

1 5+” + 

X; Gua05 







2570- 3700 

11 

1.1 

0205 



NGC 4968 

13 07 06.0 

-23 40 37 

2957 

Sy 2 

10.2 

2300- 3660 

30 

1.1 

1004 

> 10 

B; B99 

NGC 5005®) 

13 10 56.2 

-h37 03 33 

946 

LINER 

10.2 

400- 1550 

10 

1.1 

0205 

> 10 

B; Ris99 











0.3l«;= 

X; Gua05 

Mrk668 

14 07 00.4 

-h28 27 15 

22957 

RL 


22340-23630 

9 

1.1 

1004 

> 100 

X; Gua04 

Arp 220®> 

15 34 57.1 

-h23 30 11 

5434 

ULIRG 

12.1 

5250- 5720 

13 

1.1 

0601 

0.3 

C; Cle02 

NGC 6552 

18 00 07.3 

-h66 36 54 

7942 

Sy 2 

11.0 

7400- 8500 

9 

1.1 

1004 

6 + 1 

A; B99 







7400- 8500 

8 

1.1 

0205 



TXS 1946+708 

19 45 53.5 

-h70 55 49 

30228 

RL 


29000-31650 

5 

4.2 

1004 

> 28 

B; Ris03 

IRAS 20210+1121 

20 23 25.4 

-hll 31 35 

16905 

Sy 2 

11.7 

16650-17200 

8 

1.1 

0105 

< 0.6 

B; Delu04 

NGC 7674 

23 27 56.7 

-h08 46 45 

8671 

Sy 2 

11.5 

8100- 9380 

19 

1.1 

0105 

> 100 

B; B99 


a) LINER = Low Ionization Nuclear Emission Line Region; RL = Radio Loud; Sy = Seyfert; ULIRG = UltraLuminous InfraRed Galaxy 

b) Infrared flux densities taken from the IRAS Point Source Catalog (Fullmer S-c Lonsdale 1989) or from NED were converted into far infrared 
luminosities following the procedure outlined by Wouterloot Sz Walmsley (1986); Ho—75kms~^ Mpc“^; the values given are the logarithms of the 
infrared luminosity in solar units 

c) Month (first two digits) and year (last two digits) 

d) Notes concerning X-ray absorbing hydrogen column densities: X-ray telescope: A: ASCA', B: BeppoSAX; C: Chandra', X: XMM-Newton\ 
References: Arm04: Armus et al. (2004); B99: Bassani et al. (1999); Bra04: Braito et al. (2004); Cle02: Clements et al. (2002); Delu04: Deluit 
(2004); Gua04: Guainazzi et al. (2004a); Gua05: Guainazzi et al. (2005a); MaiOS: Maiolino et al. (2003); Ris99: Risaliti et al. (1999); RisOO: Risaliti 
et al. (2000); Ris03: Risaliti et al. (2003) 

e) NGC2992: a revived AGN according to Gilli et al. (2000); NGC4939: a source varying from Compton-thick to Compton-thin, see Gua05; 
NGC5005: a possibly missclassified Compton-thin source from Gua05; Arp 220: the column density from Cle02 refers to the western nucleus and is 
highly uncertain. 


overall luminosity should be fairly well correlated with 
star formation rates and thus with infrared luminosity. 
These luminosities are quite similar (^IO^^Lq) for most 
ULIRGs. Arp 220, which has slightly less than half the red- 
shift of UGC5101, should thus show ~5 times stronger 
H 2 O emission. However, Arp 220 was not detected (see 
Table 1), thus hinting at a nuclear origin of the maser in 
UGC5101. 


It is known that UGC5101 hosts a prominent buried 
AGN. Its LINER nucleus contains an absorbed hard X- 
ray component with a 6.4 keV spectral feature (Imanishi 
et al. 2003). Its central radio continuum is dominated by 
an AGN, not by a starburst (Lonsdale et al. 2003), and 
infrared signatures also point towards a dominant nuclear 
engine that is highly obscured (Imanishi et al. 2001; Armus 
et al. 2004). In the other known ULIRG hosting a water 
vapor maser, NGC6240, the H 2 O emission was found to 
spatially coincide with one of the two nuclei of the inter¬ 
acting pair of galaxies (Hagiwara et al. 2003). We thus 
conclude that the maser in UGC5101 must be nuclear. 


3.2. NGC2273 


NGC2273 is a barred, early-type spiral galaxy with a 
Seyfert 2 nucleus (e.g. Petitpas & Wilson 2002) and 
was among the galaxies observed in the Braatz et al. 
(1996) survey. The sensitivity of their measurements was, 
however, not high enough to detect a feature at the 
15mJy level. NGC2273 is by far the least luminous of 
the three detected maser sources presented here. With 
Ah 2 O< 1 OL 0 , it is a ‘kilomaser’ and may be associated 
with a prominent site of star formation instead of an 
AGN (e.g. Tarchi et al. 2002b; Henkel et al. 2005b). 
The GO emission in the innermost few arcsec of the 
galaxy (Petitpas & Wilson 2002), the nuclear dust ring 
of size 14" (Yankulova 1999), and the optical emission 
lines (e.g. Gu et al. 2003) are all indicative of ongoing 
star formation that might trigger maser emission at the 
observed luminosity. On the other hand, molecular gas 
is apparently streaming toward the innermost regions of 
the galaxy, which hosts a Compton-thick central X-ray 
source (see Sect. 5). Thus a nuclear maser is also possible. 
The linewidth, AU:^0.3kms“^ (Table 2), suggests that 
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Fig. 2. 22 GHz H 2 O maser profile smoothed to a channel 
spacing of 4.2kms“^ toward NGC2273. No emission was 
seen between 1860 and 2400 km s“^ (see Table m- 02000 
= 06*' 50“ 08!7, ^2000 = 60° 50' 45". Fsys = cz,ys = 
1840kms-i (NED). FLSR-f^HEL = +0.46kms-h 


the maser is longlived. Thus interferometric measurements 
may soon provide a position with subarcsec accuracy. 



Velocity (km/s) 


Fig. 3. 22 GHz H 2 O megamaser profile smoothed to a 
channel spacing of 16.8 km s“^ toward NGG3393. 02000 
= 10^ 48“ 23!4, 62000 = -25° 09' 43". Vtys = 3750km s"^ 
(NED). Elsr-Ehel = -9.09kms-i. 


3.3. NGC3393 

Like NGG2273, NGC3393 is a barred, early-type spiral 
with a Seyfert 2 nucleus (e.g. Gooke et al. 2000). It is 
seen nearly face-on and was also part of the Braatz et 
al. (1996) survey, where it was undetected (their Icr noise 
level is llmJy with a channel spacing of 0.66kms“^). 
More recently, however, NGC 3393 was independently de¬ 


tected by Kondratko et al. (2006). Fig.^lshows a lineshape 
reminiscent of that seen toward NGG 4258. Thus the nu¬ 
clear origin of the maser is evident, indicating the presence 
of an accretion disk that is viewed approximately edge- 
on and suggesting a rotation velocity of order 460kms“^. 
It is noteworthy that none of the disk maser (candidate) 
sources (for NGG 1068, NGG 2960, NGG 3079, IG2560, 
NGC 4945, Circinus and NGC 6323, see Gallimore et al. 
2001; Henkel et al. 2002; Kondratko et al. 2005; Ishihara 
et al. 2001; Greenhill et al. 1997, 2003; Braatz et al. 2004) 
shows the high rotation velocity, ~1000kms“^, observed 
towards NGC 4258, the prototypical source of its class. In 
view of the Toomre Q parameter for disk stability (e.g. 
Kondratko et al. 2005) it seems that NGC 4258 is hosting, 
with its comparatively large nuclear mass (e.g. Herrnstein 
et al. 1999) and circumnuclear rotation velocity (Nakai 
et al. 1993), an exceptionally compact and well ordered 
masing disk that is not typical for Seyfert galaxies. 


4. Radio versus X-ray data 

4.1. General aspects 

Traditionally, detection rates of megamaser surveys are 
low (e.g. Henkel et al. 2005b). The likely most direct way 
to identify potential megamaser sources involves a de¬ 
termination of column densities towards the nuclear en- 
gine(s). Braatz et al. (1997; their Fig. 35) proposed that 
megamaser galaxies may be sources with particularly large 
nuclear column densities. With only a dozen megamaser 
sources known at that time and with only five having 
measured X-ray absorbing column densities, however, a 
detailed study was not possible. 

While the presence of only one detected maser line does 
not strongly confine pumping scenarios, radiative trans¬ 
fer calculations provide signihcant constraints to kinetic 
temperature, H 2 density and H 2 O column density (e.g. 
Kylafis & Norman 1991). To connect H 2 O with total col¬ 
umn densities along a given line-of-sight, however, frac¬ 
tional H 2 O abundances are also needed. Such abundances 
were recently determined, not for distant AGN, but for 
some nearby interstellar clouds in the galactic disk (e.g. 
Melnick et al. 2000; Snell et al. 2000a, b, c; Neufeld et 
al. 2000a, b). The studies indicate [H 2 0]/[H2] ~ 10“®- 
10“^ in cool molecular cloud cores, ~10“® in stellar out¬ 
flow sources and ^10“"^ in warm shock heated gas. H 2 O 
column densities of order 10^® - 10®® cm“® needed for ef¬ 
ficient maser amplification (e.g. Kylafis & Norman 1991; 
Neufeld et al. 1994; Kartje et al. 1999) then require total 
hydrogen column densities well in excess of 10®® cm“®. 

With respect to observations of extragalactic H 2 O 
masers, there is considerable progress since Braatz et al. 
(1997). The sample of known extragalactic H 2 O masers 
has been quadrupled (see Henkel et al. 2005b; Kondratko 
et al. 2006). To determine column densities toward a nu¬ 
clear engine requires modeling its spectrum at X-ray wave¬ 
lengths. Seyfert 2 and LINER galaxies show nuclear X-ray 
sources with a sometimes heavily absorbed soft part of 



















































































































J. S. Zhang: Masers and column densities 


5 


Table 2. H 2 O line parameters obtained from Gaussian 
fits 


Source 

fSdV 

(Jykms”’’) 

y 

(km s~ 

AVi/2 

■') 

NGC 2273 

0.49±0.09 

1813±3 

34±9 

UGC5101 

2.69±0.13 

11721±7 

267±15 

NGC 3393 

3.96±0.66 

3151±11 

141±32 


0.75±0.20 

3720±3 

20±5 


0.71±0.22 

4250±4 

23±7 


1.70±0.22 

4309±1 

24±3 


their spectrum. Recent advances in X-ray astronomy with 
respect to angular resolution, in order to avoid contamina¬ 
tion by stellar sources and diffuse emission in the vicinity 
of the nuclear engine, provide sensible information on col¬ 
umn densities toward a significant number of AGN. 

Of particular interest is the group of ‘Compton-thick’ 
galaxies with column densities in excess of cm“^. The 
bulk of this column density must arise from a circumnu- 
clear shell or torus, since large scale (^200 pc) column den¬ 
sities in galaxies with ‘normal’ nuclear starbursts tend to 
be one to two orders of magnitude smaller (e.g. Harrison 
et al. 1999; Mao et al. 2000; Wang et al. 2004), while 
even in ULIRGs the corresponding values are <10^^cm“^ 
(e.g. Pappa et al. 2000). A large column density, of order 
10 ^"^ cm“^ or more in a highly confined region, presumably 
encompassing a few pc or less (e.g. Miyoshi et al. 1995; 
Gallimore et al. 2001), suggests that at least some fraction 
of the circumnuclear environment is dense and predomi¬ 
nantly molecular. Thus a statistical comparison between 
the occurrence and strength of nuclear maser emission and 
column density is mandatory. 

4.2. The X-ray data 

With the 53 targets from Henkel et al. (2005b; their 
Table 4), the type 2 Quasar J0804-I-3607 (Barvainis & 
Antonucci 2005), the three detections presented in Sect. 3 
and seven additional sources from Kondratko et al. (2006), 
we analyze a total of 64 H 2 O maser galaxies beyond the 
Magellanic Glouds. In order to study the correlation be¬ 
tween column density and maser emission, we searched for 
published column densities in the literature, mainly based 
on ASCA., BeppoSAX, XMM-Newton and Chandra X-ray 
spectroscopic observations. The data, including RXTE 
(Rossi X-ray Timing Explorer) observations, are summa¬ 
rized in Tables 0-0 

The sample of 64 objects includes 13 kilomasers 
(Ah 2 O< 1 OL 0 , their source names are marked by italics 
in Table 0 and 51 megamasers, the latter presumably en¬ 
tirely of nuclear origin. Except NGC4293, we find absorb¬ 
ing columns for all the kilomasers in the literature. Beside 
of NGG5194 (M51) and NGC2273, all are Gompton- 
thin. NGC5194 appears to be variable from marginally 
Gompton-thick in 1993 (Terashima et al. 1998; Bassani 
et al. 1999) to Gompton-thick in 2000 (Fukazawa et al. 
2001) and NGC2273, the newly detected maser source, is 


Table 3. X-ray flux limits for H 2 O maser sources without 
good spectra, observed and corrected for galactic absorp¬ 
tion (see Sect. 4.2). With the exception of the kilomaser 
galaxy NGC4293, all sources host megamasers. 


Source Observed Corrected 

upper upper 

limit limit 

(IQ-13 ej-gcm“^ 


NGC 235A'‘> 

2.91 

4.50 

NGC449*’> 

3.07 

6.55 

NGC591*’> 

3.73 

7.63 

NGC 613°-’ 

2.14 

3.38 

IC 184“> 

2.51 

4.19 

UGC3255'“> 

3.66 

10.34 

VII ZW 73”'> 

4.02 

9.35 

Mrk 78°-' 

5.88 

11.55 

40804+3607°' 

2.06 

3.19 

NGC 2824°' 

3.02 

5.50 

NGC 2960“' 

3.41 

6.51 

NGC 2979“' 

2.93 

5.75 

Mrk 34*’' 

1.60 

2.16 

NGC 3735“' 

1.64 

2.52 

NGC 4293“' 

2.74 

4.78 

NGC 5495“' 

4.58 

9.22 

NGC5793*’' 

9.41 

22.8 

NGC 6323“' 

1.61 

2.67 

IRAS F19370-0131“' 

5.11 

15.00 

NGC 6926“' 

3.20 

7.49 

AM2158-380NED02°' 

2.77 

4.20 

TXS 2226-184*’' 

4.31 

7.84 

IC 1481*’' 

3.61 

7.73 


a) Upper limits to 0.1-2.4 keV fluxes, obtained from the ROSAT X-Ray 
All-Sky Survey with HEASARC (High Energy Astrophysics Science 
Archive Research Center). 

b) Tentative or negative detections by ASCA (0.3-12 keV) from 
TARTARUS version 3.0 at formal 3.19, 9.34, 4.72, 7.36, 0.12, and 2.89 cr 
levels for NGC 449, NGC591 Mrk34, NGC 5793, TXS 2226-184 and 
IC 1481, respectively. For the way the given flux limits were derived, see 
also footnote a). 


a Compton-thick source detected by BeppoSAX (Maiolino 
et al. 1998) and XMM-Newton (Guainazzi et al. 2005a). 

Among the 51 megamaser sources, 27 have no pub¬ 
lished column densities. Sixteen of these were only ob¬ 
served by the ROSAT All-Sky Survey (see Table 0 foot¬ 
note a, for flux limits without and with correction for 
galactic absorption). Another six sources were measured 
by ROSAT and ASCA^ but not by BeppoSAX, XMM- 
Newton or Chandra (for limits, see Table 0 footnote b). 
Future observations of these targets by XMM-Newton and 
Chandra would be desirable. 

The remaining five sources without derived column 
densities (and not observed by ASCA or BeppoSAX) 
were recently measured by Chandra and/or XMM- 
Newton (IRAS F01063-8034, NGC 2782, ESO269-G012, 
NGC 4922, and NGC 5728). There are no published col¬ 
umn densities yet. For NGC 4922, which was observed in 
November 2004, count rates are not high enough to de¬ 
rive a meaningful result. We obtained data for the other 
four sources from the XMM-Newton Science and Chandra 
Data Archives and extracted their spectra using the 
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software packages SAS v.6.1.0 {XMM-Newton Science 
Analysis Software) and Cl AO V.2.3 (C/iondra Interactive 
Analysis of Observations, http://cxc.harvard.edu/ciao/). 
The XSPEC vll.2.0 package was then used to fit the spec¬ 
tra and to obtain column densities. 

Here we present the results of the X-ray data analysis 
of NGC 2782 and NGC 5728. The other two sources will 
be discussed elsewhere (Greenhill et ah, in preparation). 

4.2.1. Data acquisition and reduction 

Chandra observed NGC 2782 on May 17, 2002^ and 
NGC 5728 on June 27, 2003^ with effective exposure times 
of 29.6 ksec and 18.7ksec, respectively. Both observations 
were carried out with the ACIS-S detector with a stan¬ 
dard frame time of 3.2 s. The appropriate calibration 
database was provided by the Chandra X-ray Center and 
standard screening criteria were used. We checked that 
pile-up is not important due to comparatively low count 
rates of 0.03 counts/frame and 0.06 counts/frame, respec¬ 
tively, and that no flaring-background events occurred 
during the observations. Extended emission on scales of 
several arcseconds was found in both sources. Therefore, 
we extracted source spectra from circles with a radius of 
10'^ around the nuclear positions and background spectra 
from nearby regions free of X-ray sources. Nuclear spec¬ 
tra from smaller circles of 1 ” in radius were also extracted. 
Appropriate RMF- and ARF-files were created for each in¬ 
dividual spectrum by applying standard procedures. The 
energy range 0.3-8.0 keV was considered. The spectra were 
binned to contain 25 counts per bin (15 counts per bin for 
the nuclear spectra). 

4.2.2. Spatial analysis 

In Fig. 21 we show adaptively smoothed, three-color 
Chandra images of NGC 2782 and NGC 5728. Both sources 
exhibit a bright, unresolved core and extended regions 
with emission several arcseconds in size. Both large-scale 
emitting regions have optical and radio counterparts. For 
NGC 2782, the relatively weak large-scale emitting region 
is more elongated to the north than to the south, consis¬ 
tent with results from ROSAT (Jogee et al. 1998). The 
teardrop-shaped X-ray emitting region around the nu¬ 
cleus of NGC 2782 resembles the Ha structure revealed by 
Hubble Space Telescope (HST) imaging (Fig. 4 in Jogee et 
al. 1998). The Ha + Nil structure reveals the X-ray ‘knot’ 
6'.'5 south of the nucleus to be fully surrounded by a shell 
of line emitting clouds (Jogee et al. 1999). This shell is as¬ 
sociated with a starburst-driven outflow bubble. In addi¬ 
tion to a 5 GHz radio counterpart of this southern bubble, 
there is a second bubble of similar size on the opposite 
side of the nucleus, only visible at radio wavelengths. The 
detection of a hidden AGN-like nucleus from the X-ray 
spectrum of NGC 2782 (likely revealing a 6.4 keV iron fea- 

^ Obs. ID: 3014; PI: Ian Stevens 

^ Obs. ID: 4077; PI: Steven Kraemer 


ture in its innermost region; see Sect. 4.2.3) and the H 2 O 
maser emission (Braatz et al. 2004) suggest a more domi¬ 
nant role of the central engine than previously anticipated. 

NGC 5728 reveals a structure elongated southeast to 
northwest with predominantly soft X-ray emission from 
the southeast and harder emission from the northwest. 
This difference may be caused by absorption of soft X- 
rays towards the northwestern region, possibly due to a 
geometry in which the southeastern region is closer to the 
observer. This scenario is supported by optical HST emis¬ 
sion line studies that reveal a pronounced biconical ax- 
isymmetric structure of the ionization region (Wilson et 
al. 1993). The X-ray morphology resembles in size and ex¬ 
tent the inner parts of the ionization cones in NGC 5728. 

4.2.3. Spectral analysis 

Starting with the ‘canonical’ Seyfert 2 spectrum (absorbed 
power law plus soft excess; see e.g. Moran et al. 2001), 
we tried to reproduce the spectral shape of both sources 
with models of increasing complexity. It was not possible 
to achieve an acceptable fit with a simple absorbed two- 
power-law model in either case^. Not only were the re¬ 
duced chi-squared values unacceptable (3.6 and 2.4, with 
33 and 48 degrees of freedom, respectively), but the re¬ 
sulting best-fitting photon indices became unphysically 
extreme (-2.5 and 4.6 for NGC 2782, 0.6 and 9.5 for 
NGC 5728, respectively). 

We considered two alternative physical models to ac¬ 
count for the hard nuclear X-ray emission: (i) a Compton- 
thick high-energy reflection spectrum and (ii) a partially 
absorbed power law spectrum. Both models invoke a dense 
absorber towards the central source but differ physically, 
e.g. in the column density of the absorbing material. While 
model (i) implies reflection of seed photons off dense 
matter in excess of 10^"* cm“^ in a disk-like geometry, 
model (ii) applies to a scenario in which a fraction of 
the source is absorbed by intervening material (e.g., asso¬ 
ciated with dense clouds) of typically 10^^ cm“^ to a few 
times 10^^ cm“^. Model (i) has the additional benefit that 
it is sensitive to the orientation of the disk. In addition, 
we requested that a valid spectral model should fit the 
integrated source spectrum including both the hard com¬ 
pact nuclear spectrum and the spatially more extended 
soft emission. 

We modeled the soft excess via a thermal mekal 
plasma model (Mewe et al. 1985). This choice appears 
natural and physically motivated in the light of the spa¬ 
tial distribution of the soft emission (see Sect. 4.2.2 and 
Fig.lH). The transition from the integrated source spec¬ 
trum to the nuclear spectrum should thus be possible via 
modifications of the adjustable mekal parameters induc¬ 
ing variations in intensity and possibly in temperature. 


® In all our models we consider two globally absorbing com¬ 
ponents, one being representative of material inside the Galaxy, 
the other one for gas inside the measured target source. 
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Fig. 4. Adaptively smoothed three-color coded Chandra images of NGC2782 (left) and NGC5728 (right). The max¬ 
imal smoothing scales are 3" and 2", respectively. Red represents soft emission between 0.3 keV and 1.5 keV, green 
corresponds to 1.5-2.5keV and blue denotes the energy range 2.5-8.0keV. 


Table^shows the results of spectral fitting. While sta¬ 
tistically acceptable fits to the spectra of both sources can 
be obtained with both physical models, we prefer model 
(i) in the case of NGG 2782 and model (ii) in the case of 
NGG5728 (see the discussion below). Spectra and fits are 
shown in Fig.O 

The photon index of NGG 2782 in model (ii) is ex¬ 
traordinarily flat (r = 0.25). This index could only be 
explained with a complex structure of the absorbing sys¬ 
tem, invoking a variation of column densities in individ¬ 
ual absorbing clouds. While such a configuration cannot 
be excluded, the comparatively simpler model (i) yields a 
satisfactory fit to the data with reasonable values for the 
crucial parameters. The large inclination of the reflecting 
disk is consistent with a Type 2 AGN but we note that this 
parameter is, in spite of the small formal Itr errors given in 
Table 01 not well constrained by the available data. While 
the integrated spectrum does not require the addition of 
a narrow Gaussian to represent iron-line emission, excess 
emission is seen in the nuclear spectrum (radius: 1”). The 
equivalent width of this tentative feature is ^990 eV. 

In the case of NGG 5728, model (i) is statistically in¬ 
ferior to our preferred model (ii). Moreover, the high- 
energy-reflection model (i) implies a photon index as in¬ 
verted as -0.8. The partial-obscuration model (ii), on the 
other hand, yields a less extreme (although still compara¬ 
tively flat) photon index of 0.8. The large equivalent width 
of the Fe line (1.5 keV and 1.1 keV for the two models, re¬ 
spectively) is notable. 

In spite of the discussed uncertainties, the anal¬ 
ysis of the Chandra spectral data of NGG 2782 and 


NGG 5728 clearly demonstrates the presence of com¬ 
pact high-column-density absorbers in both systems. In 
NGG 2782, the high-energy-reflection model suggests a 
Gompton-thick absorber with a column density in excess 
of 10^"* cm“^ (even the alternative model (ii) still implies 
an extreme value of ~ 1 x 10^^ cm“^). For NGG 5728 the 
partial-obscuration model predicts ~ 8 x 10^^ cm“^. 

To summarize, we obtained X-ray column densities for 
38 maser sources, 36 from the literature and 2 from pro¬ 
cessing archived data. 12 of the 13 kilomaser galaxies and 
26 of the 45 megamaser galaxies are included. Because of 
the larger sample, the statistics are improved over the last 
such studies, those by Braatz et al. (1997; see Sect. 4.1) 
and Madejski et al. (2006; their Table 3 containing 11 
disk masers). 

5. Discussion 

Not all of the 38 maser sources with well defined X- 
ray spectra permit a meaningful comparison between 
H 2 O maser properties and column densities. A correlation 
would be most convincing if there were a strong dominant 
X-ray source known to coincide both with the nuclear en¬ 
gine of its parent galaxy and the observed H 2 O maser. 

These optimal conditions cannot be fulfilled. X-ray 
images with the highest angular resolution, ~1", have a 
linear resolution of ~10pc for the nearest sources (e.g. 
NGG253). For the most distant X-ray detected target, 
3G403, this angular resolution corresponds to Ikpc and 
it is hard to imagine that X-ray emission from the (hope¬ 
fully dominant) nuclear source is not contaminated by a 





J. S. Zhang: Masers and column densities 



0.5 


Channel Energy (keV) 


NGC2782 (Model I) 



Channel Energy (keV) 



Channel Energy (keV) 


NGC5728 (Model II) 



Channel Energy (keV) 


Fig. 5. Chandra X-ray spectra and spectral models of NGC2782 (left panels) and NGC5728 (right panels). The 
source spectra are from circles with a radius of 10'' around the nuclear positions. Crosses mark the observed fluxes and 
the histograms show the best-fit model, a Compton-thick high-energy reflection model for NGC2782 and a partially 
absorbed power law model for NGC 5728 (upper half of the upper panels). Residuals are given below. The lower panels 
show the unfolded model spectra with their subcomponents. In the case of NGC 2782, these are a thermal plasma 
emission (MEKAL) and a Compton-thick reflection component, modified by absorption from the Galaxy and the 
redshifted target. In the case of NGC 5728, the fit consists of a thermal plasma (MEKAL), a partially absorbed power 
law component and a Gaussian shaped Fe line, modified by the absorption from the Galaxy and the target. For details, 
see Sect. 4.2.3. 


large number of (hopefully weaker) stellar objects and an 
extended diffuse gas component. 

The maser location poses a similar problem. Towards 
some sources, the H 2 O line is observed with submilliarc- 
second accuracy. Even in the best of all cases, however, 
when H 2 O is tracing an edge-on accretion disk, the situ¬ 
ation is complex. Only the ‘systemic’ H 2 O velocity com¬ 
ponents (for a spectrum showing redshifted, blueshifted 


and ‘systemic’ components, see Fig. 0 may arise from in 
front of the nuclear engine. The red- and blueshifted fea¬ 
tures originate instead from locations displaced by up to 
a few parsec (e.g. Gallimore et al. 2001; Kondratko et 
al. 2005) from the nucleus. Amplification of the nuclear 
radio continuum by unsaturated systemic maser compo¬ 
nents is another unknown parameter that may strongly 
affect maser luminosities that depend sensitively on inch- 
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Table 4. Spectral fitting results (see Sect. 4.2.3 and Fig. El) 


Source®'^ 


Nh °'> 

kT"*) 




ph) 

< 


EW'‘> 

FluxO 

x'^/dof 




[keV] 




[keV] 

[ev] 

[keV] 



NGC2782i„t 

I 

> 10 

u.uu_o 03 

Q 1 r+O.Ol 
^•-*-^-0.01 

0.004 

75.5111 

1 Cfi+0-09 

i.t)8_Q Qg 




4.39 

52/44 


II 

12 . 7 +_ti 

^•'^'^-0.03 

^•-^'^-0.01 

0.130 


0 25 + 0-08 

- 

- 

- 

4.80 

49/45 

NGC2782nuc 

I 

> 10 

0 45+°-^^ 

1.00 

0.000 

70 4+^'^ 

2 14+0-^^ 

6 39+°'^° 


0.99 

1.86 

15/10 


II 

4 3+25 

n 90-I-0.04 

0 47+0-21 

-0.20 

0.666 

— 

~ 0.07 

6 TfS + O' 11 
D.OO_o 09 

<160 

1.32 

2.02 

15/10 

NGC5728i„t 

I 

> 10 

-0.03 

o-o27/:"i 

0.037 

^•^-2.6 

1 

0 

bo 

1 + 

0 c 

b c 

-J 0 

6-33711^ 

1207R1 

1.55 

13.0 

45/26 


II 

- 7-1-0.55 
‘ ■ ' -0.51 

0 49+°-°'^ 
^•^=^-0.03 

o.03l“;“( 

0.033 


0 77+u-*i''’ 

' ' -0.06 

fy oq + 0'33 
D.OO_o 33 

98f33 

1.13 

13.2 

35/27 

NGC5728nuc 

I 

> 10 

^•U-‘-_0.08 

^•'^^-0.02 

0.005 

Q 7+34.9 
' -9.7 

±.ut?_o 07 

6 ‘12+0-0'^ 
0-'5-^_0,03 

10873? 

1.29 

9.54 

23/15 


II 

°-^-0.50 

n 1 r-l-0.01 
^•-*-^-0.01 

0 14+0.03 
^•-^‘*-0.03 

0.011 


0 94+0-08 
'-'•^^-0.07 

6 ‘12 + 0-04 
'^•'^^-0.04 

851^9 

0.89 

9.34 

16/16 


a) sourceint: Integrated spectra from circular areas with a radius of 10^^ around the nuclear position; sourceNuc: Nuclear spectra from smaller areas 
of l'' in radius 

b) M: Model I: a Compton-thick high-energy reflection spectrum (wabs*zwabs(mekal-|-hrefl(pow-Kgau)) in XSPEC); Model II: a partially absorbed 
power law spectrum (wabs*zwabs(mekal-hpow-hzwabs(pow-|-gau)) in XSPEC). The galactic absorbing column densities are 1.76X 10^°cm~^ and 
7.83x lO^^cm”^ for NGC2782 and NGC5728, respectively (from Dickey & Lockman 1990) 

c) X-ray absorbing column density in units of lO^^cm”^ 

d) Temperature of the thermal MEKAL plasma in keV 

e) Metallicity of the thermal MEKAL plasma as a fraction of solar metallicity (in the case of model I of NGC2782 nuc 5 this relative metallicity is 
not constrained by the data and was therefore fixed to unity). 

f) Ratio between the directly seen continuum and the scattered/absorbed continuum 

g) Disk inclination 

h) Photon index 

i) There is no obvious line in the integrated spectrum of NGG 2782 (NGC2782int)- Intrinsically present line emission may be buried in the relatively 
bright extended continuum 

j) Width of the line feature that is tentative in the case of NGG 2782 (its strength does not significantly affect the other parameters) 

k) Equivalent width of the line feature, see footnote j 

l) Flux from 0.3 to 8 keV in units of 10~^^ergcm~^s~^ (corrected for galactic absorption) 


nation and warping of an approximately edge-on maser 
disk. 

There are further complications. For many of the 
sources with high column density, only lower limits of 
order cm“^ have been derived (see TableEJ. 

There is a major uncertainty in the column density esti¬ 
mates related to the possibility that the X-ray emission 
is predominantly scattered light not passing through the 
main body of the obscuring torus. A significant fraction 
of the kilomasers are related to star forming regions well 
off the center of their parent galaxy. In IC 342, neither the 
maser nor the dominant X-ray source are nuclear (Tarchi 
et al. 2002a; Bauer & Brandt 2004), both being associated 
with different star forming regions. In NGC 598 (M 33) the 
‘nuclear’ X-ray continuum originates from an X-ray binary 
close to the dynamical center, well off the H 2 O masers 
that are observed in the large-scale disk. There is even 
one H 2 O detected galaxy, M82, that lacks a well defined 
nuclear engine. 

All these complications suggest that a correlation be¬ 
tween H 2 O maser emission and X-ray column density 
should not only show a large scatter, but may remain 
entirely undetectable. For the 38 maser sources with 
well defined X-ray spectra, column densities (ignoring all 
the problems outlined above, preferring high resolution 
Chandra data whenever possible and adopting the most 


recent value if Chandra data are not available) are dis¬ 
played in Fig.EJi.. The histograms show the number of ex- 
tragalactic H 2 O masers as a function of column density, 
both for the entire sample and for the subsamples of 12 
kilomasers and 26 megamasers. We find that there is a 
striking difference between the column density distribu¬ 
tion of the two classes, with megamasers being associated 
with higher column densities. 

All kilomasers (see TableEl their source names are 
printed in italics) except NGC 2273 and NGC 5194 are 
Compton-thin, i.e. their absorbing column densities are 
NH<10^‘‘cm“^. There are two groups of sources, one with 
column densities of 10^^“^^ cm“^ and another one with 
column densities >10^^'^ cm“^. The mean of the logarith¬ 
mic values is <log A^h (cm“^)> = 22.62 ± 0.37, the error 
being the standard deviation of the mean (for NGC 2273, 
the lower limit to the column density was taken; see 
TableEJ . 

H 2 O megamasers exist preferentially in environments 
with higher obscuration, i.e. in circumnuclear regions of 
high column density. The peak in the distribution lies at 
24 <log Vh (cm“^)< 25. 85% (22/26) of the megamaser 
sources are heavily obscured (VH>10^^cm“^) and one half 
are Compton-thick (A^H^lO^'^cm”^). These fractions are 
consistent with the average properties of Seyfert 2 galax¬ 
ies as derived by Risaliti et al. (1999) and Guainazzi et 
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al. (2005b), while showing slightly but not significantly 
higher column densities than those reported by Bassani 
et al. (1999). Taking logarithmic values, the mean col¬ 
umn density becomes <logA^H (cm“^)> = 23.77 ± 0.14. 
For six sources, there are only lower limits; these adopted 
limits thus lead to an average below the true mean 
value, likely underestimating the difference between mega¬ 
maser and kilomaser column densities (see TableEJ. The 
Kolmogorov-Smirnov statistical method was used to test 
whether the distributions of column density of our kilo¬ 
maser and megamaser sources are statistically distinguish¬ 
able (see Fig.Eb). The K-S test result shows that the null 
hypothesis, they both are from the same parent popula¬ 
tion, can be rejected at a confidence level of at least 99.4%. 

In view of the uncertainties outlined above and the fact 
that, to date, no systematic searches for H 2 O in galaxies 
with Compton-thick nuclear environments have been un¬ 
dertaken, the result is surprisingly clear: Nuclear column 
density is a very solid guide when searching for mega¬ 
maser emission. Almost all megamasers known to date 
could have been found by analyzing X-ray data. 

The cause of the difference between kilomaser and 
megamaser column densities is less clear. Differences in 
position between maser and X-ray sources were already 
mentioned. On the one hand, we may face a high degree of 
small scale clumping in the observed clouds, even within 
circumnuclear accretion disks. On the other hand, mil- 
liarcsecond offsets between megamaser and nuclear X-ray 
source are expected, while such offsets tend to be much 
larger in some of the kilomaser sources. Even accounting 
for the smaller mean distance to known kilomaser galax¬ 
ies, linear offsets between maser and X-ray sources tend 
to be much larger in kilomaser sources than for mega¬ 
masers. The three kilomaser sources, however, that are 
known to be ‘nuclear’ (i.e. the maser and dynamical cen¬ 
ters coincide within 1") and possess an associated X-ray 
source (NGC253, NGC4051, and NGC5194) reveal col¬ 
umn densities that are indistinguishable from those of the 
megamaser sample. While three sources are clearly not 
sufficient for a statistically meaningful result, we note that 
this effect alone has the potential to explain the discrep¬ 
ancy in averaged column densities between kilomaser and 
megamaser column densities. 

The subsample of (candidate) disk masers with X-ray 
column densities, encompassing NGG1068, NGC2639, 
NGC3079, IC2560, NGG 3393, NGG 4258, NGG 4945, 
and Circinus, is another subsample worth to be studied. 
For the mean of the logarithmic values of their column 
densities we find <log A^h (cm“^)> = 23.96 ± 0.36. This 
value (two sources are represented by lower limits only) 
is slightly higher but still consistent with that derived for 
the entire megamaser sample. 

In Fig.UI the H 2 O maser isotropic luminosities are 
plotted against their X-ray absorbing column densities. 
The H 2 O maser isotropic luminosities were taken from 
Table 4 of Henkel et al. (2005b) and from Sect. 3. We 
adopted a ±30% uncertainty in maser luminosity for all 
sources. For the entire maser sample, again taking lower 


limits to the column densities as real values for seven tar¬ 
gets, the linear fit gives logLnaO = (0.56±0.17) log A^h 
± (-11.89T3.93), with a correlation coefficient r=0.49. A 
trend of rising column density with maser luminosity may 
be apparent, but the scatter is large and the correlation 
is weak. The megamaser subsample (not displayed) shows 
no trend at all. There is no exponential growth of inten¬ 
sity with column density as expected in the case of un¬ 
saturated maser emission (e.g., Goldreich et al. 1972). A 
correlation roughly fitting TH 2 O oc A^^ would be expected 
in the case of idealized saturated maser emission (no ve¬ 
locity gradients in the maser region) with the value of the 
exponent being determined by the luminosity increasing 
linearly with column density and the surface of the masing 
cone growing with the square of its lengths (e.g., Kylafis 
& Norman 1991). As already indicated, even slight differ¬ 
ences in the lines of sight as well as amplification by a 
background continuum source may be responsible for the 
weakness of the correlation. 

6. Conclusions 

In this paper, maser emission and column density are dis¬ 
cussed for 64 extragalactic H 2 O maser sources. The main 
results are: 

(1) In searching for the A~1.3cm (22 GHz) water va¬ 
por line in galaxies with high nuclear column densities, ac¬ 
tive nuclei or high infrared luminosities, we detected H 2 O 
emission towards NGG 2273, UGC5101, and NGG 3393 
with isotropic luminosities of 7, 1500, and 400 Lq, respec¬ 
tively. UGC5101, one of the most luminous H 2 O mega¬ 
masers ever found, is the first luminous (Lh 2 O> 1 OOL 0 ) 
H 2 O maser detected in an ULIRG. The maser is likely 
nuclear. NGG 2273 hosts a maser that is either nuclear 
or associated with massive star formation. Since its line- 
shape is wide, suggesting a long lifetime, interferomet¬ 
ric measurements may soon reveal its nature. NGG 3393, 
first detected by Kondratko et al. (2006), shows a clas¬ 
sical NGG 4258-like profile. If the analogy holds, the ro¬ 
tation velocity of the disk is well below the 1000 km s“^ 
encountered in NGG 4258. This and similarly low rotation 
velocities in other ‘disk-maser’ galaxies suggest that the 
accretion disk in NGG 4258, surrounding a nuclear engine 
of relatively high mass, is particularly compact and well 
ordered. 

(2) Chandra X-ray images and spectra are presented 
and analyzed for NGG 2782 and NGG 5728. A hid¬ 
den AGN was detected in NGG 2782. The X-ray spec¬ 
trum is best fit by a high energy reflection model with 
Ain^lO^^cm”^; its 6.4keV Fe line appears to be buried 
in a spatially more extended, relatively bright continuum. 
The H 2 O maser, so far believed to be associated with star 
formation, may be related to the nuclear engine. The X- 
ray spatial analysis of NGG 5728 confirms the biconical 
axisymmetric structure of its ionization region, which was 
previously revealed by optical HST emission line studies. 
The best fitting model, with a partially absorbed power 
law spectrum, yields A^h 8 x 10 ^^ cm 
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Fig. 6. a. Histograms showing the number of H 2 O masers 
(kilomasers: Lh 2 O< 1 OL 0 ; megamasers: Lh 2 O> 1 OL 0 ) as 
a function of column density determined from X-ray spec¬ 
troscopy (upper panels). Because of angular resolution, 
Chandra data are taken whenever possible. Otherwise 
the most recent measurements are used (see Table 0 ). 
b. Kolmogorov-Smirnov test for the kilomaser and mega¬ 
maser samples (lower panel). The maximal difference of 
the cumulative probability (the vertical line) is 0.564 and 
the probability that one parent population would produce 
the kilomaser and megamaser galaxy samples is 0.006. 


(3) A statistical analysis of the correlation between X- 
ray absorbing column density and H 2 O emission was car¬ 
ried out, encompassing 64 galaxies beyond the Magellanic 
Clouds. At a confidence level of at least 99.4%, the col¬ 
umn density distributions of kilomaser and megamaser 
sources are different. More than 80% (22/26) of the mega¬ 
maser sources are heavily obscured (A^H> 10 ^^cm“^) and 
half are Compton-thick (A^h >10^'^cm“^). Among the 13 
known kilomaser sources (TH 2 O < 10 L 0 ), however, only 
NGC2273 and NGC5194 are Compton-thick. The small 
subsamples of nuclear kilomasers and disk megamasers 
show average column densities that are undistinguishable 
from those of the entire megamaser sample and from sam¬ 
ples of Seyfert 2 galaxies that were not selected with 
respect to maser emission. Presumably a clumpy cloud 
structure in the cirumnuclear environment, diverging posi- 



LogN^ 


Fig. 7. Isotropic maser luminosities (in L 0 ) versus X-ray 
absorbing column densities A^h (in cm“^), including a lin¬ 
ear fit. The data are from Table 0 taking Chandra data 
whenever possible and adopting otherwise the most recent 
result. The six arrows in the upper right corner represent 
seven sources: Mrk 1066 and NGC 5347 show similar H 2 O 
luminosities (Henkel et al. 2005b) and lower limits to the 
X-ray column density (Table0. 


tions between maser and nuclear sources and an occasional 
amplification of a background radio continuum source are 
sufficiently decoupling X-ray column densities and H 2 O 
maser properties not to show a clear correlation. 

(4) The identification of an AGN of type 2 is sufficient 
to justify a search for H 2 O maser emission. The determina¬ 
tion of X-ray column densities toward the nuclear engine 
does not enhance maser detection probabilities. 
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Table 5. Nuclear X-ray absorbing column densities of extragalactic H 2 O maser galaxies* 


Source 

Type'^ 

Telescope^ 

Epoch® 

Nh 

References*^ 





( 10 ^^ cm“^) 


ICIO 

X-1 

c 

2003.03.12 

'^•'J'J-0.0008 

BB04 

NGC235A 

S2 

R 


see Table 3 

HEASARC 

NGC253 

SBC 

B 

1996.11.29 


Cap99b 



C 

1999.12.16 


Wea02 

NGC262 (Mrk348) 

S2 

A 

1995 


Ris02 

IRAS F01063-8034®) 

S2 

X 

2003.04.21 

see Sect. 4.2 


NGC449 (Mrkl) 

S2 

A 

1996.08.06 

see Table 3 

TARTARUS 

NGC598 (MSS) 


A 

1993.07 

0 - 016 t°;oo 2 

Tak94 


X -8 

X 

2000-2002 

0.019t°;°!]i 

Fos04 

NGC591 

S2 

A 

1997.07.28 

see Table 3 

TARTARUS 

NGC613 

S 

R 


see Table 3 

HEASARC 

IG184 

S2 

R 


see Table 3 

HEASARC 

NGC1052 

LINER 

A 

1996.08 

0 

^•'^-0.83 

Ter02 



X 

2001.08.15 

~ 7 

Kad06 

NGC 1068 

S2, C 

A 

1993.07.24 

> 100 

B99, Iwa97 



B 

1996, 1998 

> 100 

Mat97, Gua99 



C 

2000 . 02.21 

> 100 

Mat04 



X 

2000.07.29 

> 100 

Mat04 

Mrk 1066 

S2, C 

A 

1997.07.24 

> 10 

Ris99 

NGC 1386 

S2, C 

B 

1996.12.10 

> 10 

Mai98 



X 

2002.12.29 

> 22 

Gua05 

IG 342 


A 

2000.02.24 

0.084tO;°« 

Kub02 


X-21 

X 

2001 . 02.11 

0.087t°;°i3 

Kon03 

UGC3255 

S2 

R 


see Table 3 

HEASARC 

Mrk 3 

S2, C 

A 

1993 

13lf 

Tur97b 



B 

1997.04.16 


Cap99a, MattOO 



X 

2000.10.19 

13.6t°i 

Bia05 

NGG2146 

SBC 

A 

1997.03.26 

u.uzi_Q gii 

Del99 

VIIZW 73 

S2 

R 


see Table 3 

HEASARC 

NGG2273 

S2, C 

B 

1997.02.12 

> 100 

Mai98 



X 

2003.09.05 

> 18 

Gua05 

Mrk 78 

S2 

R 


see Table 3 

HEASARC 

Mrk 1210 

S2, C? 

A 

1995 

> 10 

B99 



B 

2001.05.05 

-1 qq+0.08 

Ohno04 



X 

2001.05.05 

2 14 + 0.20 

Gua02b 

J0804-b3607 

Type 2 Quasar 

R 


see Table 3 

HEASARC 

NGC 2639 

LINER 

A 

1997.04.16 

4.2+2'3 

B99, Wil98 

NGC 2782®) 

SBG, C 

C 

2002.05.17 

> 10 

this paper 

NGC 2824 (Mrk 394) 

S2 

R 


see Table 3 

HEASARC 

NGC 2960 (Mrk 1419) 

LINER 

R 


see Table 3 

HEASARC 

UGC5101 

ULIRG 

C 

2001.05.28 

0.1 (<0.65) 

Ptak03 



C,X 

2001 

131" 

Arm04 

NGC 2979 

S2 

R 


see Table 3 

HEASARC 

NGG3034 (M82)^^ 

SBG 

B 

1997.12.06 

0.058t°;°i4 

Cap99b 



X 

2001.05.06 

0.017+O;Oi« 

Ste03 



C 

2002 

0.012 

Str04 

NGC 3079 

S2/LINER, C? 

A 

1993.05.09 

0 17+018 
n-i'- 0.12 

B99, Ter02 



B 

2000.05.26 

~ 100 

lyoOl 



C 

2001.03.07 

n 17+0.02 
n-i'- 0.02 

Cec02 

Mrk 34 

S2 

A 

1997.11.17 

see Table 3 

TARTARUS 

IC2560 

S2 

A 

1996.12.19 

~3 

IshOl 



C 

2000.10.29 

> 10 

Iwa02 



C 

2004.02.16 

^30 

Mad06 
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Table 5. continued. 


Source 

Type'^ 

Telescope^ 

Epoch® 

Nh 

References*^ 





(10^^ cm“^) 


NGC3393 

S2, G 

X 

2003.01.05 

44^23 

Gua05 

NGC3556 

X-35 

c 

2001.09.08 

0 n‘?+ooi7 

Wang03 

Arp 299 (NGC3690) 

SBG, G 

X 

2001.05.06 

~ 26 

Ball04 



G 

2001.07.13 


Zez03 



B 

2001.12.14 

or 0+1.39 
^^•^-0.56 

Del02 

NGG3735 

S2 

R 


see Table 3 

HEASARG 

NGC4051 

Sl/1.5 

A 

1993.04.25 

U.UUO_Q Q35 

Geo98 



G 

2000.04.25 

~ 0.01 

GolOl 



X 

2001.05.16 

~ 3.6 

Pou04 

NGG4151 

S1.5 

A 

1993.12.07 

U.DOO_q 035 

Geo98 



A,B 

1999-2000 

2 .01°;} 

Sch02a 

NGG4258 

S1.9 

A 

1999.05.15 

0.95l°;2i 

ReyOO 



B 

1998.12.19 

1 oi +0.07 
-‘-•^-‘--0.06 

Ris02 



G 

2000-2001 

~ 0.7 

You04 



X 

2000.12.08 

q q 
d d 
+ 1 
O 
00 

o 

Pie02 



X,G 

2000-2002 

~ 1.0 

Fru05 

NGG4293 

LINER 

R 


see Table 3 

HEASARG 

NGG4388 

S2 

A 

1995.07.21 

3.34lJ;°9° 

For99 



B 

2000.01.03 

4.80ll«° 

Ris02 



X 

2002 .12.12 

' ^-0.07 

Bec04 



G 

2001.06.08 

q r-l-0.4 

I-wa03 

ESO269-G012®) 

S2 

G 

2004.01.05 

see Sect. 4.2 


NGG 49228) 

S2 

G 

2004.11.02 

see Sect. 4.2 


NGG4945 

S2, G 

A 

1993.08.31 

40l};° 

B99, Don96 



B 

1999.07.01 

22 l° 

GuaOO 

NGG5194 (M51) 

S2, G 

A 

1993.05.11 


B99,Ter98 



B 

2000.01.18 

selfe 

FukOl 

NGG5256 (Mrk266) 

S2 

A 

1999.05 

0 .161°;?° 

LevOl 

NGG 5347 

S2, G 

A 

1997 

> 10 

Ris99 

NGG 5495 

S2 

R 


see Table 3 

HEASARG 

Gircinus 

S2, G 

A 

1995.08.14 

> 10 

B99, Mat96 



B 

1998.03.13 

43l? 

Mat99 



G 

2000.06.06 

~ 60 

SamOl 

NGG 5506 (Mrkl376) 

SI? 

B 

1998.01.14 

A Q7+0.05 
-0.05 

Ris02 

NGG 5643 

S2, G 

B 

1997.03.01 

> 100 

Mai98, B99 



X 

2003.02.08 

6-10 

Gua04 

NGG 5728®) 

S2 

G 

2003.05.27 

7 7+0.55 
'’^-0.51 

this paper 

NGG 5793 

S2 

A 

1998.08.05 

see Table 3 

TARTARUS 

NGG 6240 

ULIRG, G 

A 

1994.03.27 

>20 

I-wa98 



B 

1998.08.15 

21 .81?-? 

Vig99 



X 

2000.09.22 

lo.oi?;? 

Bol03 



X 

2001.07.29 

13.01^0° 

Ptak03 

NGG 6300 

S2, G? 

RX 

1997.02 

> 100 

Lei99 



B 

1999.08 

2 in+°-°° 

Gua02a 



X 

2001.03.02 

9 1 c-l-O.OS 

Mats04 

NGG 6323 

S2 

R 


see Table 3 

HEASARG 

ESO 103-G035 

S2 

A 

1996.03.18 

2 lf)+0-26 

For99 



B 

1997.10.14 

o ao+0.28 

Z.UZ_o 28 

WilOl 

IRAS F19370-0131 

S2 

R 


see Table 3 

HEASARG 

3G403 

FRII 

G 

2002.12.07 

4 00“'"°-^° 

U.UU_Q 20 

Kra05 

NGG 6926 

S2 

R 


see Table 3 

HEASARG 

AM2158-380NED02 

S2 

R 


see Table 3 

HEASARG 

TXS 2226-184 

LINER 

A 

1997.11.21 

see Table 3 

TARTARUS 
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Table 5. continued. 


Source 

Type'^ 

Telescope^ 

Epoch‘s 

iVn 

(10^^ cm“^) 

References*^ 

IC 1481 

LINER 

A 

1997.11.28 

see Table 3 

TARTARUS 


* Chandra and XMM-Newton results are given whenever possible. For sources not observed by these satellites most 
recent (but usually older) results are quoted. Sources in italics: kilomaser sources. 

a) Type of nuclear activity. In the absence of a starburst and a prominent nuclear source, the X-ray source (e.g. 
X-1 or X-8) is given. SBG: StarBurst Galaxy; S2: Seyfert 2; LINER: Low-Ionization Nuclear Emission Line Region; 
ULIRG: UltraLuminous InfraRed Galaxy; ERIE Eanarov-Riley Type II galaxy; G: Gompton-thick, i.e. 

C?: possibly varying between Gompton-thick and Compton-thin 

b) A: ASCA; B: BeppoSax; C: Chandra] R: ROSAT] RX: RXTE] X: XMM-Newton 

c) Year, month and day 

d) Arm04: Armus et al. (2004); Ball04: Ballo et al. (2004); B99: Bassani et al. (1999); BB04: Bauer & Brandt 
(2004); Bec04: Beckmann et al. (2004); BiaOS: Bianchi et al. (2003); Bia05: Bianchi et al. (2005); Bol03: Boiler et al. 
(2003); Cap99a: Cappi et al. (1999); Cap99b: Cappi et al. (1999); Cec02: Cecil et al. (2002); ColOl: Collinge et al. 
(2001); Del99: Della Ceca et al. (1999); Del02: Della Ceca et al. (2002); Don96: Done et al. (1996); For99: Forster et 
al. (1999); Fos04: Foschini et al. (2004); Fru05: Fruscione et al. (2005) FukOl: Fukazawa et al. (2001); Geo98: George 
et al. (1998); Gua99: Guainazzi et al. (1999); GuaOO: Guainazzi et al. (2000); Gua02a: Guainazzi (2002); Gua02b: 
Guainazzi et al. (2002); Gua04: Guainazzi et al. (2004b); Gua05: Guainazzi et al. (2005a); IshOl: Ishihara et al. (2001); 
Iwa97: Iwasawa et al. (1997); Iwa98: Iwasawa & Comastri (1998); Iwa02: Iwasawa et al. (2002); Iwa03: Iwasawa et 
al. (2003); lyoOl: lyomoto et al. (2001); Kad06: Kadler et al. (2006), in prep.; Kon03: Kong (2003); Kra05: Kraft et 
al. (2005); Kub02: Kubota et al. (2002) Lei99: Leighly et al. (1999); LevOl: Levenson et al. (2001); Mad06: Madejski 
et al. (2006); Mai98: Maiolino et al. (1998); Mat96: Matt et al. (1996); Mat97: Matt et al. (1997); Mat99: Matt et 
al. (1999); MatOO: Matt et al. (2000); Mat04: Matt et al. (2004); Mats04: Matsumoto et al. (2004); Mil90: Miller & 
Goodrich (1990); Ohno04: Ohno et al. (2004); Per98: Persic et al. (1998); Pie02: Pietsch et al. (2002); Pou04: Pounds 
et al. (2004); Ptak03: Ptak et al. (2003); ReyOO: Reynolds et al. (2000); Ris99: Risaliti et al. (1999); Ris02: Risaliti 
et al. (2002); SamOl: Sambruna et al. (2001); Sch02a: Schurch et al. (2002a); Sch02b: Schurch et al. (2002b); Smi96: 
Smith & Done (1996); Ste03: Stevens et al. (2003); Str04: Strickland et al. (2004); Tak94: Takano et al. (1994) Ter98: 
Terashima et al. (1998); Ter02: Terashima et al. (2002); Tur97a: Turner et al. (1997a); Tur97b: Turner et al. (1997b); 
Vig99: Vignati et al. (1999); Wang03: Wang et al. (2003); Wea02: Weaver et al. (2002); WilOl: Wilkes et al. (2001); 
Wil98: Wilson et al. (1998); You04: Young & Wilson (2004); Zez03: Zezas et al. (2003). 

e) Absorbing column densities were obtained by processing XMM-Newton and Chandra Archive data (see Sect. 4.2). 
}) For NGC3034 (M82), lacking a well defined nucleus, the Chandra data refer to the diffuse halo component. 

g) Chandra archive data have not yet been released. 



